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Abstract 

We predict a clean signal at the Large Hadron Collider (y / s=14TeV) 
for a scenario where there is a top- like, charge +| vectorlike isosin- 
glet fermion. Such a quark, via mixing with the standard model top, 
can undergo decays via both flavour-changing Z-boson coupling and 
flavour-changing Yukawa interactions. We concentrate on the latter 
channel, and study the situation where, following its pair-production, 
the heavy quark pair gives rise to two tops and two Higgs boson. We 
show that the case where each Higgs decays in the bb channel, there 
can be a rather distinct and background-free signal that can unveil 
the existence of the vectorlike isosinglet quark of this kind. 

1 Introduction 

Our present knowledge about elementary particles and their interactions upto 
the energy scale of several hundred GeV's is encapsulated in the theory called 
the Standard Model(SM). The SM, a renormalizable gauge theory of strong 
and electroweak interactions based on the guage group SU(3)c x SU(2)l x 
U(1)y, gives a successful explanation for most of the phenomena governing 
fundamental processes, and is in excellent agreement with the experimental 
data to date. However, there are a number of unanswered questions which 
motivates us to think beyond the SM. These include, just to name a few, the 
flavour and naturalness problems, the absence of a cold dark matter candi- 
date in the spectrum, and the origin of neutrino masses and mixing pQ. They 
have led to a plethora of conjectures extending the SM. 
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On the experimental side, a golden opportunity to test many of these conjec- 
tures has come through the Large Hadron Collider (LHC), which aims at not 
only discovering the only missing piece of the SM, namely the Higgs boson, 
but also looking for new physics beyond it. 

One of the various ideas beyond the SM is to extend the fermionic sector of 
the SM by postulating existence of either a sequential fourth generation [2] or 
vector-like fermions. The electroweak precision data strongly constrains the 
existence of extra chiral fermions. On the other hand, vector-like fermions, 
have left-and right-handed components with the same transformation prop- 
erty under SU(2), and are considerably free from the aforementioned con- 
straints. They can have gauge invariant mass terms of the form i/jl^r, which 
do not arise from the electroweak symmetry breaking mechanism, and should 
be traced to some new physics scale. Thus their signature is of immediate 
interest if that scale is accessible to the LHC, with the fond hope that they 
might contain some clue to the flavour problem. 

Vectorlike fermions appear as singlets, doublets or triplets under SU(2), in 
many options beyond the Standard Model(BSM) like Little Higgs Model[3], 
composite Higgs model [1] and extra dimensional models [5]. They also appear 
in some grand unified theories like -E^IE], which once got impetus from consid- 
erations underlying superstring theories. In particular, vectorlike isosinglets 
also play a role in the AF = 2 effective Lagrangian analysis studies[7]. Also, 
it has been pointed out that models with such an extended quark sector 
can give rise to quark electric dipole moment at the two-loop level due to 
the lack of Glashow-Illiopoulos-Maiani(GIM) suppression, thereby implying 
constraints on the model parameter space(s)[8]. The collider phenomenology 
of such an isosinglet vector fermion has been studied from various angles 

Here we focus on an SU(2) singlet, charge | vectorlike quark. It should be 
noted that down-type vectorlike isosinglets, too, have been considered ex- 
tensively in the literature [8], [29]- [32]. The vectorlike quark is pair-produced 
at the LHC in the same way as the top quark, subject, of course, to the 
inevitable kinematical suppression if its mass is higher. Such a quark, how- 
ever, has additional decay channels, which can make it distinct. The main 
features responsible for such distinction are its isosinglet character, and its 
capacity to mix with the top, once SU(2) x U(l) is broken. As a result of 
doublet-singlet mixing in the left chiral sector, such a quark, named t' here, 
has flavor changing interactions both with the higgs boson(H) and the Z 
boson. We propose to utilise the resulting dacay channels, namely t'— > tH 
and tZ. In particular, we find that the former channel leads to a rather 
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unique and background free signature arising from a ttHH state. 

On account of mixing between the SM fermions and their SU(2) singlet coun- 
terparts, many observables are expected to be different from the Standard 
Model predictions, specially in the sector involving third family. We assume 
a situation where Higgs is discovered already and its mass is approximately 
known. We use such a Higgs as an instrument to investigate the production 
and decay of exotic top-like quark into channels which hardly have Standard 
Model backgrounds. We consider two Higgs masses, namely m# = 120 GeV 
and ran — 130 GeV. In order to maximize our signal rates, we let the Higgs 
decay in each case into the final state with maximum branching ratio, namely 
H — > bb. We show that, using the consequent 6b final states(including de- 
cays of the top quarks as well), one can construct signals with very little SM 
backgrounds. Tagging five b's out of six in each case, with appropriate event 
selection criteria, proves to be sufficient for this purpose. We demonstrate 
that such a signal can allow us to probe in a discriminating fashion a large 
part of the parameter space consisting of the t'mass and its mixing with the 
t. The rest of the paper is organised as follows. Section 2 contains an outline 
of the scenario, statements on the signal looked for, a reminder of the con- 
straints on various parameters, and a resume of the methodology adopted. 
The results are presented in Section. 3. We summarise and conclude in 
Section 4. 

2 The scenario, the signal and the methodol- 
ogy 

2.1 The scenario and its signals 

As has been said already, we consider a minimal extension of the SM with 
the inclusion of a top-ike vector isosing let, t' L (3,1,+f), t' R (3,1,+f) to the 
matter content of the SM. The gauge boson and the Higgs sector remains 
unchanged. 

t' can be produced in pair via the strong interactions or singly via eletroweak 
processes. We concentrate here on the former channel, which at the parton 
level corresponds to gg — > t't'and qq — > t't! essentially arising from the gluon 
coupling of the heavy quark: 

Lqcd = -ig s Vi^ G » C 1 ) 

Neglecting the small contribution from electroweak diagrams, the pair pro- 
duction cross section for an isosinglet fermion and a chiral fourth generation 
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fermion is the same [53] . [34] . The production cross section depends only upon 
the mass of t'and goes down with increase in the mass. 

Once the 577(2) x U(l) symmetry is broken, the most substantial mixing of 
t'can take place with the top, as there are rather stringent bounds on mix- 
ing with the first two generations |14j. It must be remarked that there have 
been studies which considered the mixing with the lighter generations, which 
can affect, for example, the single production of vector like quarks through 
the electroweak channels |35j. As a result of such mixing, the four charge | 
quarks in the weak basis are related to the corresponding mass eigenstates 
by 
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V is the Standard Model CKM matrix. In such a basis the mass matrix 
for the down type quarks is diagonal. The mass matrix(M u )for the up type 
quarks is 



L lL<lR 
M f T n 



F t' 



(3) 



where q^R = (u,c,t)L,R, Mg LqR is 3 x 3 mass matrix of the SM particles, 
Mn T fi is 3 x 1, Mn, t i is 1 x 3 and Mp *> is the mass term for t'. Mj „ 
and Mp p do not arise from the Yukawa couplings. M u is diagonalized by 

L ft 

the bi-unitary transformation : WM U U' = M% iag , where U' is parametrized 
analogously to U but in spite of right handed fermions being there in the 
Yukawa couplings, the elements of U' do not appear as such. 
With the structure of the mixing matrix being what is shown in equation(T5]), 
the Standard Model CKM matrix is no longer unitary, and instead forms a 
block in the unitary 4x4 mixing matrix U. The V and W together form 
the 4x3 charged current mixing matrix. The violation of CKM unitarity 
leads to a breakdown of the Glashow-Illiopolous-Miami(GIM) mechanism, 
and leads to the flavor changing neutral current (FCNC) processes (t —> cZ) 
and (t — > cH) in the top sector [101 [TTJ. As mentioned above, t' decays to 
the SM fermions along with either the electroweak guage bosons (W ,Z) or 
Higgs boson(H) at the tree level. The charged current interaction (£'-> bW + ) 
is given by 

L cc = ^ff L ^b L W; (4) 
On account of mixing with the top quark we have Flavor changing neutral 
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current interactions with Z boson given by 



J neutral 



2Cos 



t' L l^ L Z^ (5) 



And the interaction of t'with Higgs boson and the Flavor changing Yukawa 
coupling is 

Lyukawa = -yt>qLiH c t' R + h.c (6) 

where y t > in this comes to be V^Wm M# . In addition to the above yukawa 
coupling there are terms proportional to t' L qn and t' L t' R which arise on account 
of the fact that it is not possible to diagonalize the mass and Yukawa matrices 
simultaneously. 

Following the simplified version of the mixing matrix used in [11] , we describe 
all the interactions of t' by the following mixing matrix 



U 
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(7) 



Where its interactions with all quarks of the first two generations are ne- 
glected. We also neglect the interactions of the SM quarks across the gen- 
erations. The elements of the mixing matrix are chosen with the motivation 
to have considerable flavor changing decay modes in addition to the charge 
changing decay modes. 

We concentrate on one of the decay process of t', namely, (£' — > tH). For 
the sake of our current work we assume that the Higgs boson has already 
been discovered with approximately known mass. This particular decay 
channel has already been used as a tool to study the discovery potential 
of Higgs [T3l [261 [27] but with a different final state. In [28] though the final 
decay of h considered is also to bb but it is in context of "Little Higgs model" 
where there are extra gauge and Higgs bosons along with t' . The further 
decay chain we consider is where both the tops decay to bW and both the 
Higgs to bb. i.e 

pp t'i' hthi -»■ bbbW + bbbW~ (8) 

The final state consists of 6b's and 2 Ws, out of which we attempt to identify 
5b's and predict the signal for 5b+X final state. As we show in this work, 
this turns out to be a rather clean signal with very low SM background. 

Direct and indirect searches for the vector fermions put constraints on the 
mass and couplings of t'and the mixing angle 9, between t and t' . Among 
them, direct searches imply the following bounds: 
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• Considering t' — > W + q as the only possible decay chain, the lower limit 
on the mass of t'is 358 GeV at 95% CL. by the CDF collaboration [3"3~| 
at center of mass energy, a/s=1.96 TeV 

• The DO collaboration puts the lower limit in the channel of W+jets 
decay, at the 95% CL to be 258 GeVjS] at ^=1.96 TeV. 

• The latest study by ATLAS set the lower limit on the mass of t'at 
m t i < 404 GeV at the 95% CL. assuming 100% decay through bW + 
mode|36] at the center of mass energy, \fs=l TeV. 

• Assuming 100 % branching fraction for the decay if— >tZ, t'with any 
mass less than 475 GeV is excluded at 95 % confidence level by CMS 
detector at the LHC{37] at ^/s=7 TeV. 

• The R b ratio given by, R B = r{ z^drlns) > § ives a stron S constraint on 
the t-f mixing angle, 9 to be < 25° [IE]. 

For the purpose of our current work, we consider the lowest mass to be 350 
GeV, and make no specific assumption about branching ratios in the three 
decay channels. 



2.2 Methodology 

For numerical calculations, we have used CalcHEP v2.5.6[3H] and a CalcHEP- 
PYTHIA interface program^ along with PYTHIA-6.4.240Q]. The pro- 
duction cross-section of the isosinglet quark pair and the decay of if , if is 
calculated using CalcHEP. In order to do so, a new model with the new in- 
teractions based on the mixing matrix considered, was added to the existing 
list of CalcHEP models. We have taken m t =172 GeV, and used CTEQ6L 
parton distribution functions (PDF), for the center-of-mass energy a/s=14 
TeV. We chose the following benchmark points for our calculation: After 



Parameter 


Value 


M v [Gev\ 


350,400,500 


mixing angle t — t', 6 


5,10,15 


M h [GeV] 


120,130 



Table 1: Benchmark Points 



calculating the production cross section and branching fractions for all the 
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benchmark points by CalcHEP, the output is transferred to the CalcHEP- 
PYTHIA interface program [39j. The interface program is provided with the 
appropriate selection cuts. Subsequently, PYTHIA is used for obtaining rates 
in the desired final state of 5b+X. 

The SM background in this case is calculated at the parton level for the signal 
pp — > hhtt using CalcHEP v2.5.6(3B]- The only on shell process which gives 
rise to this signal in the SM is pp — > ti — > hhtt, as a result the background 
is very small, b-identification efficiency has been taken to be 50%. We have 
not taken into account the effect of, for example, charm-induced jets faking 
b's. We reconstruct the invariant mass of all the bb pairs which survive the 
first three of the selection cuts listed below, and follow the criteria explained 
there, the background for the signal 5b+X is very low in comparison to the 
signal. 

We have used the following selection criteria on the minimum of five b's 
required in our stipulated final state: 

• Each of the identified b's should have E T >40.0 GeV. 

• Each b jet should be central, with pseudorapidity, | r\ |<2.5 . 

• We implement b-tagging efficiency of 50% i.e e& ~0.5 . 

• As a final step we calculate the invariant mass for all the possible 
combinations of b pairs. We impose the following restriction on the 
calculated invariant mass(mbf,) i.e m&& = {Mu ± 15) GeV. 

We get our final numbers by counting all the events(NH) which have 
at least two such b pairs with their calculated invariant mass falling in 
the above limit. We predict our signal using this number. 

3 Results 

We are presenting the results at the leading order (LO), and thus our esti- 
mate can be called conservative. Our final state consists of six b's and two 
Ws. A similar use of the final state comprising six b's has been considered 
in [20]. However, the suggested signal, clean as it is, has suppressed rates 
due to the requirement of ^identification of all the six b's and ii) the simul- 
taneous tagging of an isolated lepton. We, in contrast, give up the lepton 
tagging requirement. Moreover, we suggest identifying only five out of six 
b's, with the proviso that four out of them display two individual peaks, 
each at the mass of the Higgs boson which is presumably detected before our 
analysis takes place. We succeed in considerably enhancing the event rates 
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in this manner, while at the same time having negligible backgrounds, with 
our chosen acceptance criteria. As is clear from the plot, figure (TJ the pair 
production cross section of t' decreases with M t > since it depends only on the 
mass. We find that the production of signal at the parton level, hhtt in the 
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Figure 1: Production Cross section of t!t! with the mass of t' at \/s=14 TeV 

SM, is weaker by atleast 10 2 order of magnitude in comparison to the lowest 
signal(for M t '=500 GeV, M^=130 GeV and 9 = 15) in our model as can be 
seen in the table ([2]) and (EJ). We present our results as the cut-flow chart 



M h (GeV) 


Xsection[pb] 


120 


0.00060 


130 


0.00044 



Table 2: Production Cross section in the Standard Model for the signal 
pp hhtt at y/s=14 TeV 





M h [GeV] 


Cross-section 


pb] 


9 = 5 


9=10 


9 = 15 


350 


120 


0.4272 


0.4105 


0.3826 


130 


0.3590 


0.3583 


0.3355 


400 


120 


0.2430 


0.2365 


0.2232 


130 


0.2308 


0.2236 


0.2101 


500 


120 


0.07089 


0.06931 


0.06520 


130 


0.06949 


0.06761 


0.06324 



Table 3: Production cross section of pp — > tthh in our model at -y/i=14 TeV 

for all the considered benchmark points and the invariant mass distributions 
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M t ,[GeV] 


Cut 


Cross-section [pb] 


M h =120 GeV 


M h =130 GeV 


6 = 5 


9=10 


6=15 


6 = 5 


6=10 


6=15 


350 


E b T > AO.OGev 


0.1993 


0.1907 


0.1809 


0.1730 


0.1664 


0.1575 




Vb 


< 2.5 


0.1735 


0.1660 


0.1574 


0.1519 


0.1462 


0.1378 


NH > 2 


.01061 


.01046 


.0098 


.0099 


.0094 


.0089 


400 


E b T > AO.OGev 


0.1155 


0.1120 


0.1061 


0.1072 


0.1038 


.0973 




Vb 


< 2.5 


0.1023 


.0951 


.0939 


.0954 


.0925 


.0868 


NH > 2 


.00606 


.00593 


.00557 


.0060 


.0060 


.0054 


500 


E b T > AO.OGev 


.0338 


.0323 


.0304 


.0321 


.0311 


.0292 




Vb 


< 2.5 


.0307 


.0293 


.0275 


.0293 


.0283 


.0267 


NH > 2 


.00175 


.00165 


.00155 


.0017 


.0016 


.0016 



Table 4: Outflow table for various benchmark points for M^=120 GeV and 
M h =130 GeV at ^=14 TeV 



of bb pairs. We compute our results for different values of m' t and Mh and 
mixing angle 6. In spite of a substantial reduction of the signal due to tag- 
ging efficiency of 50% per b, our predicted signal is still good enough to be 
observed at the LHC at the integrated luminosity of 30 fb" 1 . We find the 
following trends in our results as can be seen from the table HI 

• For a given value of M[ and Mh change in 6 does not make much of a 
difference( .2fb-.6fb) [figure H E|. 

• For a given value of Mh and 6 dependence on the M t ' is the strongest. 
It changes the cross section for the signal by 8fb-10fb. 

• For a given M[ and 6 the change in Higgs mass does not really make a 
difference. 

• Also as the M' t goes from 350 GeV to 500 GeV, the change in Mh and 
6 hardly makes any difference on the signal cross-section. 

We plot the invariant mass distribution of all the combinations of b pairs for 
the events which survive the first three selection criteria. We find that this 
distribution has two peaks corresponding to the reconstructed Higgs mass 
with in the required mass range (m&6 = {Mh ± 15)GeV), as in figure !4][ol 
satisfying our fourth selection criterion, for all our benchmark points. We 
present here the distribution for two points only i.e. fig SI [5]). 
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0.01 - 

0.009 - 

0.008 - 
% 0.007 - 
m 0.006 - 

a 0-005 - 

0.004 - 

0.003 - 

0.002 - 

0.001 -- 
30 

Figure 2: Cross section for the final state (5b+X) with the mass of t'for 
various values of mixing angle 9 and Mh—120 GeV at - v /s=14 TeV 

M h =130 GeV 

0.01 -j 
0.009 - 
0.008 - 
0.007 - 
* 0.006 - 
0.005 - 
0.004 - 
0.003 - 
0.002 - 
0.001 -- 

300 350 400 450 500 

M,[GeV] 

Figure 3: Cross section for the final state (5b+X) with the mass of t'for 
various values of mixing angle 9 and M^=130 GeV at -^5=14 TeV 

4 Summary and Conclusions 

We have calculated the particular signal (5b+X), in the framework of a model 
with a top like(+|) vector fermion, [if), in addition to the SM particles. We 
assume that it mixes only with the top quark. We ignore not only the 
interactions of t'with first two generations but also the interactions of the 
SM quarks across the generations in the mixing matrix, U in equd As a 
result of the mixing with the top, t' has Flavor Changing Neutral interactions 
with Z and H bosons. We make use the interaction with Higgs boson and 
follow a particular decay mode (t' — > ht) of t'and further consider the decay 
of higgs to hi. Using CalcHEP, PYTHIA and CalcHEP-PYTHIA interface 
programs, we predict an observable signal of 5b+X from a final state signal 
of 6b's and 2Ws, at the LHC at the center of mass energy, - v /s=14TeV. We 



M h =120 GeV 




350 400 450 500 

M,[GeV] 
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Figure 4: Invariant mass distribution of bb for r?ty=350 GeV, mixing angle 
9=15 and M h =120 GeV for y/s=U TeV 



Figure 5: Invariant mass distribution of bb for m f /=500 GeV, mixing angle 
9 = 15 and M h =120 GeV for y/s=U TeV 



consider the situation when the Higgs boson is already discovered, so that 
its mass is a known quantity, and can be used to identify two b-pairs with 
invariant mass around the mass of the Higgs, taken here to be in the region 
120 GeV-130 GeV. We find that in spite of the rather ambitious proposal of 
tagging five b's, we get, after all cuts, a signal of the order of few fb which 
is stronger by at least 10 2 in comparison with the SM background. Since all 
our results are at the leading order so the predictions we make about the 
signal are rather conservative. We conclude that an integrated luminosity of 
30 fb' 1 should be sufficient to either find out or rule out the existence of £'in 
the mass range 350 GeV-500 GeV. 
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